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SLIDE  1 


IC  Is  generally  agreed  that  an  important  component  of  the  sea  surface- 
generated  radiation  is  a  diffuse  distribution  of  vertical  dipole  sources. 1,2,3 
Such  radiators  result  from  spray  Impact  on  the  surface  and  the  injection  of 
nonspherical  air  bubbles  a  small  fraction  of  a  wavelength  from  the  pressure 
release  surface.  The  dipole  source  structure  is  fundamental  to  nearly  all  of 
the  acoustic  ambient  prediction  algorithms  in  current  use.  Despite  the  ubi¬ 
quitous  use  of  the  term  there  are  relatively  few  measurements  of  the  source 
strength*'3'8  and  thus  little  is  known  about  its  statistical  properties.  It 
was  the  objective  of  the  task  being  discussed  here  to  estimate  the  statistical 
properties  of  these  dipole  source  strengthr. 

The  approach  taken  in  this  study  had  two  components.  First,  a  data  base 
was  accumulated  from  measurements  of  the  solid  angle  vertical  directional 
spectrum  made  in  an  acoustically  isolated  area  where  the  acoustic  ambient  is 
dominated  by  the  local  sea  conditions.  Second,  the  dipole  source  strength 
“area  density"  is  found  by  matching  only  that  part  of  the  total  measured 
directional  spectrum  that  exhibited  a  dipole  pattern. 
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SLIDE  2 


The  Tongue  of  the  Ocean  in  the  Bahamas  was  the  location  of  the  measure¬ 
ment.  It  is  a  relatively  deep  water  basin  totally  isolated  from  global  ship¬ 
ping.  The  actual  cul-de-sac  measurement  location  is  further  shielded  from 
transiting  vessels  north  of  the  site  by  the  landmass  north  of  the  location. 
There  is  no  place  to  transit  to  in  the  cul-de-sac  itself  so  shipping  is  nearly 
nonexistent  in  the  vicinity  of  the  sensor  system.  The  acoustic  measurement  is 
some  30  nmi  from  the  anemometer  which  is  10  m  above  the  water  surface. 

The  data  base  consists  of  101  trials  taking  place  over  a  ten-month 
period.  The  minimum  time  between  trials  is  48  hours  which  for  the  most  part 
assured  statistical  independence. 


SLIDE  3 


The  measurement  system  objective  was  a  broadband  characterization  of  the 
diffuse  acoustic  ambient  field.  This  slide  shows  the  acoustic  sensor  system 
used.7  It  consisted  of  a  vertical  hydrophone  array  which  covered  the  two 
decade  frequency,  range  of  40  to  4000  Hz  with  seven  octavely  nested  antennas 
each  with  a  modest  2  to  4  wavelength  aperture.  A  parametric  spectral  estima¬ 
tion  algorithm  that  modeled  the  media  refraction  and  source  directional  char¬ 
acteristics  was  then  used  to  obtain  higher  angular  resolution  near  the  hori¬ 
zontal  direction.7 

The  seven  octaves  were  recorded  sequentially  by  synchronously  recording 
the  eight  14-bit  hydrophone  signals  from  each  antenna  system  on  video  cassette 
recorders.  The  depth  and  tilt  of  the  vertical  hydrophone  array  were  also 
recorded. 

The  brightly  colored  syntactic  foam  collar  around  the  instrumented  pres¬ 
sure  vessel  maintains  the  verticality  of  the  system  within  a  small  fraction  of 
the  hydrophone  array  angular  resolution. 


INSTRUMENTATION  SUMMARY 
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SLIDE  4 


This  slide  summarizes  the  procedure  used  to  obtain  the  dipole  source  area 
density  function.7  The  first  step  in  the  data  processing  is  to  estimate  all 
possible  cross  spectral  density  (CSD)  functions  between  hydrophone  outputs. 
These  terms  are  the  elements  of  the  CSD  matrix  on  the  left-hand  side  of  the 
first  equation.  The  left-hand  side  of  the  equation  is  the  weighted  quadratic 
form  of  the  CSD  matrix  which  calculates  the  beamformed  output.  Equivalently, 
the  right-hand  side  is  also  the  beamformed  output  found  by  the  solid  angle 
integral  of  the  product  of  the  ambient  directional  spectrum  and  the  hydrophone 
array  pattern  function.  The  procedure  was  to  model  the  directional  spectrum 
as  a  sum  of  terms  with  unknown  linearly  occurring  parameters  which  are  then 
evaluated  so  as  to  minimize  the  mean  square  difference  between  the  left  and 
right-hand  sides  of  the  equation. 

The  form  of  the  model  is  ad  hoc,  physically  based,  and  all  of  the  compo¬ 
nents  of  the  model  occur  historically  in  the  acoustic  ambient  literature.7 
For  our  present  purpose  there  are  three  requirements  of  the  model.  First,  the 
sum  of  the  components  must  adequately  represent  the  measurements.  Second,  the 
components  must  be  distinguishable  from  each  other,  and  third,  one  of  the 
terms  must  model  a  diffuse  field  of  vertical  surface  dipoles.  The  latter  term 
is  shown  explicitly  in  the  slide.  It  consists  of  the  independent  addition  of 
all  possible  sources  contributing  to  a  specific  elevation  angle  at  the  obser¬ 
vation  point  per  unit  solid  angle.  The  radiation  from  each  source  is  propa- 
gaced  to  the  observer  with  a  quite  general  solution  of  the  wave  equation 
incorporating  all  pertinent  bottom,  surface,  and  water  column  terms.  The  par¬ 
ameter,  aD,  is  the  source  area  density  evaluated  by  linear  least  squares 
estimation  problem.  It's  dimension  is  jipa2  per  Hz  per  m2. 
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PARAMETRIC  SPECTRAL  ESTIMATION  SUMMARY 
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T( )  is  propagation  term  between  xs  and  x. 


SLIDE  5 


This  linear  axis  figure  illustrates  the  general  character  of  the  dipole 
source  frequency  spectrum.  The  abscissa  is  frequency  and  the  ordinate  is 
spectrum  level.  These  are  three  examples  of  the  over  100  such  spectra  found. 
The  spectrum  are  clearly  bandpass  with  a  decaying  high  frequency  end  and  an 
abrupt  low  frequency  cutoff.  The  shape  is  reminiscent  of  gravity  wave  spectra 
and  the  data  smoothing  functions  illustrated  are  the  same  analytical  form  used 
to  model  gravity  wave  spectra  but  with  quite  different  parameters.  The  envi¬ 
ronment  is  parametized  by  the  nondimensional  surface  friction  velocity  used  by 
B.R.  Kerman;3  i.e.,  the  surface  friction  velocity  is  nondijnensioned  by  the 
minimum  phase  speed  of  the  capillary- gravity  wave  field.  Unity  value  of  this 
coefficient  represents  the  onset  of  clearly  visible  spilling  whitecaps. 
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SLIDE  6 


Repeating  the  previous  slide,  only  this  time  with  the  more  familiar 
logarithmic  scales,  allows  an  examination  of  the  low  frequency  character  of 
the  data.  While  the  high  frequency  spectra  is  uniformly  inversely  propor¬ 
tional  to  frequency  with  an  exponent  typically  slightly  greater  than  one  the 
low  frequency  falloff  rate  varies  widely.  Sections  of  the  spectra  can  be  fit 
with  2nd  and  3rd  power  exponents.  No  regularity  of  these  low  frequency  shapes 
with  the  environmental  parameter  could  be  found. 
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SLIDE  7 


In  this  figure  example  dipole  source  area  spectrum  are  compared  with  pre¬ 
vious  measurements.11'5'6  The  log-log  axis  presentation  is  the  same  as  the 
previous  slide  but  the  environmental  parameter  is  changed  to  be  the  10  m  wind 
speed  with  examples  chosen  to  match  the  previous  data.  The  comparison  indi¬ 
cates  two  points.  First,  this  data  generally  agrees  with  previous  estimates 
when  spilling  breakers  are  numerous.  This  is  particularly  true  with  the  Kup- 
erman  and  Ferla6  data.  Second,  when  spilling  breakers  are  not  present  the 
abrupt  falloff  in  the  spectrum  occurring  in  the  present  data  is  not  observed 
in  the  previous  measurements . 
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This  figure  quantifies  the  first  order  description  of  the  bandpass  char¬ 
acter  of  the  data.  The  abscissa  is  the  logarithm  of  the  nondimens ional  fric¬ 
tion  velocity  and  the  ordinate  is  frequency  in  Hz.  The  figure  shows  the 
frequency  at  which  each  of  the  measured  spectrum  peaked  and  the  upper  and 
lower  half- power  points.  It  is  seen  that  while  the  upper  frequency  changes 
little,  both  the  frequency  of  peak  energy  and  the  lower  half-power  frequency 
are  inversely  proportional  to  the  surface  friction  velocity  with  a  break  point 
at  incipient  spilling  breakers.  Numerically  the  actual  bandwidth  is  nearly 
1400  Hz,  independent  of  the  environmental  variable. 


t 
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NONDIMENSIONAL  FRICTION  VELOCITY 
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The  final  figure  displays  the  frequency  integral  of  the  previous  spec¬ 
trum;  i.e.,  the*  total  dipole  source  level  in  /ipa2  per  m2.  This  variable  is 
displayed  on  a  logarithmic  scale  on  the  ordinate.  The  abscissa  is  again  non- 
dimensional  friction  velocity.  Also  indicated  in  the  figure  is  uB  which  is 
identified  by  E.C.  Monahan  and  M.  Lu®  as  the  lowest  values  of  the  friction 
velocity  for  which  there  is  bubble  injection  from  the  water  surface  into  the 
water  column.  This  is  a  10  m  wind  speed  value  of  about  2  to  2.5  m/s.  For 
friction  velocity  between  Uj  and  uc  the  data  is  highly  variable  but 
monotonically  increasing.  Taking  the  attitude  that  the  source  level  should 
increase  with  bubble  injection  rate,  an  analytical  form  typically  fit  to 
whitecap  coverage  data  was  used  in  this  region.  Above  uc  the  now  compact  data 
was  readily  fit  with  a  power  law.  With  this  in  mind  we  see  that  the  figure  is 
consistent  with  the  idea  that  there  are  three  regimes  present.  First,  there 
are  no  dipole  sources  present  below  uB/uc.  Second,  the  dipole  source  level  is 
proportional  to  the  bubble  injection  rate  until  the  onset  of  visible  spilling 
breakers.  The  source  level,  like  the  bubble  injection  rate,  is  highly  vari¬ 
able  when  parametized  by  the  surface  friction  velocity  indicating  that  other 
near- surface  dynamics  are  active.  In  fact,  whitecap  coverage  may  well  be  the 
better  independent  variable  in  this  range.  Once  the  nondimens ional  friction 
velocity  exceeds  unity  the  variability  decreases  dramatically  and  the  dipole 
source  level  becomes  proportional  to  the  surface  wind  stress. 
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SLIDE  10 


The  final  figure  is  a  summary. 
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Dipole  Source  Frequency  Spectrum 
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Dipole  Source  Area  Density  Function  Model  Presented 
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